Dark Energy Perturbations Revisited 
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In this paper we study the evolution of cosmological perturbations in the presence of dynamical 
dark energy, and revisit the issue of dark energy perturbations. For a generally parameterized 
equation of state (EoS) such as wn(z) = wo + wiy^ , (for a single fluid or a single scalar field 
) the dark energy perturbation diverges when its EoS crosses the cosmological constant boundary 
u>d = — 1. In this paper we present a method of treating the dark energy perturbations during 
the crossing of the wd = — 1 surface by imposing matching conditions which require the induced 
3-metric on the hypersurface of vjd = —1 and its extrinsic curvature to be continuous. These 
matching conditions have been used widely in the literature to study perturbations in various models 
of early universe physics, such as Inflation, the Pre-Big-Bang and Ekpyrotic scenarios, and bouncing 
cosmologies. In all of these cases the EoS undergoes a sudden change. Through a detailed analysis of 
the matching conditions, we show that 8d and 9d are continuous on the matching hypersurface. This 
justifies the method used [JQ in the numerical calculation and data fitting for the determination 
of cosmological parameters. We discuss the conditions under which our analysis is applicable. 

PACS numbers: 98.80.Cq; 95.36.+X 

I. INTRODUCTION 

Since the discovery that the expansion of the universe has recently been accelerating, a discovery made in particular 
through observations of distant Type la supernovae (SNIa) in 1998 @, @, a lot of effort has been made to understand 
the reason for the acceleration. The most popular interpretation of the data is to assume that the current universe 
is dominated by a new form of matter with negative equation of state denoted "dark energy" . The equation of state 
(EoS) wd of the dark energy, defined as the ratio of its pressure to energy density, is usually used to classify the 
different dark energy models. One of the candidates for dark energy is the cosmological constant whose EoS wd is 
a constant and ecjuals —1 at all times. In dynamical dark energy models extensively discussed in the literature, such 
as quintessence [JQ, phantom (Toj. k-essenceflll]. anintom(l2j and so on, wd is generally a function of the redshift 1 . 
For quintessence dark energy —1 < Wjj, while for phantom wd < — 1- The salient feature of quintom dark energy is 
that its EoS crosses the phantom boundary set by w = — 1. 

Given this wealth of theoretical models for dark energy, it is crucially important to use the accumulated high 
precision observational data from SNIa, Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) 
surveys to constrain the value of wd and its evolution. In this data-driven investigation, one needs to begin by 
parameterizing wd(z) and then fit the parameters introduced to the data. In recent studies, a popular parametrization 
of the EoS is the CPL parametrization [20l [2l| : Wd(z) = wq + w\ , where wo and w% arc two free parameters. 
This model is simple and has a clear interpretation: wq is the present value of the EoS and u>i is its derivative with 
respect to the scale factor a. In Figure 1 we show the whole parameter space of this model. Interestingly it can be 
divided into four regions by the two blue dotted lines which according to the classifications of the models in terms 
of the EoS correspond to quintessence, phantom and quintom A and quintom B respectively. Both quintom A and 
quintom B have Wd crossing — 1. However, they cross the cosmological constant boundary in a different way. For 
quintom A, Wd transits from wd > — 1 to wd < ~ 1 as the universe expands, but quintom B does in a opposite way. 
The crossing point of the two blue dotted lines corresponds to the model of the cosmological constant. 

The EoS merely reflects the nature of dark energy on the background of a homogeneous universe. Unless we restrict 
our attention to the special case of the cosmological constant, we must take into account the dark energy perturbations 
to obtain a consistent and complete procedure for data analysis, in particular when fitting cosmological parameters 



1 For example, see Refs. [13H19I ] for reviews on dark energy. 
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to the data of CMB and LSS. In fact, it has been shown that the results obtained by data fitting are quite different 
depending on whether one includes or docs not include the dark energy perturbations (for examples, see, fH [22T[2a |). 

With the energy and momentum density perturbations of dark energy denoted by Sd and 6d, the perturbation 
equations are simple if we assume that dark energy consists of a single perfect fluid or a single scalar field. In this 
case, if is restricted not to cross the line of —1, the perturbation equations behave well. However, such an a-priori 
restriction on the EoS will yield a biased result because it excludes most of the parameter space of the model. Thus, 
in order not to loose generality, one should do the global data analysis for the whole parameter space as shown in 
Figure 1. However, when wd crosses the line —1, the perturbations will diverge [l2l. [27l - l29j . In fact, in this context of 
General Relativity as the theory of space-time, it is impossible to obtain a background which crosses the "phantom 
divide" with only a single scalar field or a sin gle p erfect fluid. This is why the quintom scenario of dark energy needs 
to introduce extra degrees of freedom 0, H2I l27l43l| 2 . And this also implies that the parametrization of wd(z) for 
the background evolution will not be applicable anymore when considering the perturbations consistently. In turn, it 
makes the data analysis more complicated and inconvenient. In order to keep the maximal generality with the least 
free parameters for the parameterized EoS, Ref. [l[ proposed a method to deal with the dark energy perturbations 
during the crossing of the boundary wd ~ — 1 3 . In this method a small positive parameter e is introduced which 
divides the whole time interval into three regions corresponding to times when wd > —1 + e, when wd < —1 — e and 
when wd is between — 1 — e and — 1 + e, i.e. the region when it crosses —1. In the regions with wd > — 1 + e and 
wd < — 1 — e the perturbation equations can be solved easily. In the region when — e < 1 + wd < e, <5d and Bd are 
taken to be constant so that they are continuous in the whole time range. 

In this paper, we will revisit the issue of dark energy perturbations and will pay particular attention to the treatment 
of the perturbations when they cross the line wd = —1 taking a different view of point. Our starting point is the 
general relativistic matching conditions across space-like hypersurfaces [36l l37j which are in turn generalizations of 
the Israel matching conditions across time- like hypersurfaces |38| . These matching conditions tell us how the metric 
and its first derivative are related on the two sides of a distributional source of matter which leads to the transition 
between one solution of General Relativity on one side of the surface to a different solution on the other side of the 
surface. 

We consider the space-like hypcrsurfacc u>£>(?7, xl ) = — 1 (with r/ denoting conformal time and x l the spatial 
coordinates). The matching conditions tell us that the induced 3-metric on this hypersurface and its extrinsic curvature 
should be continuous across the matching surface [36l l37j . These matchin g co nditions have been used widely in studies 



of perturbations in various models of early universe, such as inflation [361 137| , pre-big-bang cosmology [39| , Ekpyrotic 
cosmology [4(| Sl[ and non-singular bouncing cosmologies [42l - |45j in which the EoS undergoes a sudden change. 
Through the analysis of matching conditions, we will show in this paper that - under certain conditions which will be 
discussed later - Sd and 9d are indeed continuous on the matching hypersurface, thus justifying the method suggested 

in Q- 

The present paper is organized as follows: in Section II we briefly review the difficulties encountered in single field 
or single fluid dark energy models when the EoS crosses — 1, the solution to this problem obtained by introducing 
quintom model, and the approach of Ref. [l[ on how to deal with quintom dark energy perturbations when fitting 
to observational data. In Section III we then study the transfer of dark energy perturbations across the phantom 
transition from the point of view of matching conditions on the hypersurface of wd = —1- Using the results of this 
method we then perform a global analysis of the EoS of dark energy within the framework of the CPL parametrization, 
making use of current cosmological observations. Our numerical results show that the dark energy perturbations 
cannot be neglected. In Section IV we summarize and discuss our results. 



II. GENERAL CONSIDERATION OF THE DARK ENERGY PERTURBATION 

In our analysis, we pick the Conformal Newtonian gauge in which (in the case of a spatially flat universe), the 
metric including linear fluctuations is given by 

ds 2 = a(?y) 2 [(l + 2<P)di] 2 - (1 - 2^)S l3 dx l dx j ] , (1) 

(we are focusing on scalar perturbations only 4 ). The metric perturbations $ and depend on space and time 
and describe the small deviations from a homogeneous Friedmann-Robertson- Walker (FRW) universe. They are 



2 For a consistent and complete proof of the no-go theorem, please see, [22l |. 

3 See e.g. I32I435H for relevant study on the dark energy fluctuations through the cosmological constant boundary. 

4 We refer to Ref. |4(J for a comprehensive review of cosmological perturbation theory. 
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determined by the matter perturbations through the Einstein equations which take the following form when expanded 
to linear order 

-fc 2 ^-3-H(*' + 'H$) = AnGa 2 Sp , 

k 2 (V' + H$) = 47rGa 2 (p + p)6 , 
k 2 

+ H(2^' + <P') + (2n' + -H 2 )$ + — (* = AnGa 2 Sp , 

fc 2 (*-$) = 12nGa 2 (p + p)a , (2) 

where % = din a/dry is the conformal Hubble parameter and the prime denotes the derivative with respect to conformal 
time. The energy density and pressure perturbations are denoted by Sp = STq and Sp = —1/3577, respectively. The 
variable 8 denotes the momentum density perturbation, which is defined by 

(p + p)6 = ik l ST? . (3) 

The shear perturbation a relates to the anisotropic stress through the relation 

(p + p)a = kiktfTj-l/MjSTl), (4) 

and it vanishes if matter is a perfect fluid or consists of a set of scalar fields as in the cases considered in this paper. 
Thus, in the cases considered here we have ^ = 3>. Given these considerations, one can obtain the Poisson equation 
from the Einstein equations ([2]): 

k 2 TW 

^$ = -4irGp{6 + + w)6] , (5) 

where S = Sp/ p is the density contrast. 

If there are many components of matter, then each species has its own perturbation variables Si, Oi and Spi. The 
total perturbations are given by the sum over of all species: 

pS = y]pi$i , 

i 

(p + p)9 = Y^(pi +Pi)8i , 

i 

Sp = ^2Spt . (6) 

i 

If there are no interactions beyond gravitational ones among these components, the perturbations for each species 
satisfy the individual energy and momentum conservation laws [13] 

s[ = -(i + Wl )(e t -W)~m( 6 -^~wA 

\ Pi 

e[ = -H(l 3^)0, - + k 2 ( + $) . (7) 

l + Wi \ 1 + w l J 

To solve these two equations we need to know how Spi depends on <5, and 0^: 

S Pl = c 2 J Pl + m(l + Wi)%{c 2 si - c 2 J9 t , (8) 



where 



c: 
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w l ~w'J[3H{l + w l )} (9) 



is called adiabatic sound speed in the literature, and c S i is the sound speed defined in the comoving frame of the fluid. 
For a perfect fluid c S i = c a i and for a canonical scalar field c S i = 1. 

For the problem discussed in this paper, we assume that the universe is filled with only two components, the non- 
relativistic matter (including dark matter and baryons) and the dark energy. The matter perturbation equations are 
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well behaved. But the dark energy perturbation diverges when its EoS crosses —1. For example, if the dark energy 
is a single scalar field, we have the perturbation equations from ((7|) 

8 D = -(i + WD )(6 D - W) - m{cl D WD) s D - m w ' D + m{1 ~ Wd) D , 

e' D = -n(i - 2,c 2 sD )e D + ^^5 D + fc 2 $ . (io) 

1 + WD 

When wd crosses —1 and if the sound speed remains positive, then the second equation of (fTU)) becomes singular and 
9' D diverges. For a single fluid, when crossing —1 the sound speed 

4d = c Id =wd- w' d /[3H(1 + wd)} (11) 

is divergent and can be arbitrarily negative. This is another way to see that if the dark energy consists of a single 
degree of freedom, its EoS cannot cross the boundary —1, otherwise the perturbation equations become singular and 
lead to gravitational instability. Thus, to realize the quintom scenario we should introduce extra degrees of freedom. 
The simplest quintom model is constructed by a combination of a quintessence field and a phantom field [13 ] . We 
know that the perturbations of a system including both quintessence and phantom are stable, so in the quintom model 
there is no gravitational instability. 

The quintom scenario with multi fluids or multi fields of dark energy can allow wd to consistently cross the 
cosmological constant boundary. However, it introduces more parameters for data fitting. Here, we would like to 
keep the number of free parameters for the parameterized EoS the same as for a single fluid model. A technique 
was developed in Ref. [l| to treat the perturbations during the time interval when the dark energy equation of state 
crosses the line Wd = — 1 with the goal of applying the procedure to quintom models. Specifically the authors of 
Ref. [l[ introduced a small positive parameter e. For regions with wd > — 1 + e and wd < — 1 — e, the dark energy 
behaves like quintessence and phantom, respectively, and the perturbations can easily be evolved. In the region of 
— e < 1 + wd < e, Ref. [l[ assumed 5 and 9 to be constant, i.e., 

s' D = o e' D = o (12) 

during this phase. Another way of stating this assumption is that the values of Sd and 0d are matched between the 
two sides of this interval, that is 



Sd\- d \+ = D \- , (13) 



where 5d\± and 9d\± represent the corresponding values when 1 + wd = ±e. Hence, Sd and 9d are continuous 
throughout. With this method, in fitting the data one does not need to introduce more parameters. The numerical 
calculations have shown that this approach approximated results obtained using quintom models to a high precision 
for values of the parameter e as small as 10 -5 [H-Q. 



III. DARK ENERGY PERTURBATION WITH PARAMETERIZED EOS AND MATCHING 

CONDITIONS 

The method to deal with dark energy perturbation with its EoS across —1 proposed in [![ assumes that in the 
neighborhood of the crossing point the energy and momentum density perturbations 5d and 9d are frozen. This 
guarantees the continuity of Sd and 9d- In this section we will investigate this treatment from a different point of 
view. 

Consider a space-like hypcrsurface £ which divides space-time into the two regions wd > w + = — 1 + e and 
wd < W- = — 1 — e. The surface represents the region — 1 + e > wd > — 1 — e, the region in which the evolution 
of dark energy fluctuations is not under control for EoS with a single component. To this surface we apply the 
matching conditions of [36l . which state that the induced 3-metric on this hypersurface and its extrinsic curvature 
are continuous. 

These matching conditions can be applied to reheating in inflationary cosmology: instead of solving the equations 
of motion in a specific model which describes the transition from the inflationary phase to the radiation phase after 
reheating, we cut out a time interval i_R — e < t < tji + e about the reheating time tn and apply the matching conditions 
to connect the fluctuations on either side of this interval. Similarly, these matching conditions have been applied to 
pre-big-bang and Ekpyrotic cosmology to cut out a time interval about the time when the background is singular and 
then connect the fluctuations on either side of the matching surface. 
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As pointed out in j39| . this matching procedure is not well justified if the background does not obey the matching 
conditions. Therefore, in recent studies of non-singular bouncing cosmologies [481] one introduces a bouncing phase 
valid around the bounce point and matches both at the boundary between the initial contracting phase and the onset of 
the bouncing phase, and then once again between the end of the bouncing phase and the final expanding phase. In this 
case, the fluctuations can also be evolved numerically and one can verify that the approximate analytical description 
of the evolution of fluctuations using matching conditions gives accurate results for the evolution of cosmological 
perturbation. Note that in this case the equation of state of the background also has w crossing — 1 at the bounce 
point HI. 

However, in our present investigation the matching prescription for fluctuations is justified since the background 
satisfies the corresponding conditions. Thus, it is sufficient to use a single matching surface, like in the case of 
inflationary reheating. 

In a homogeneous universe, the matching hypersurface coincides with that of fixed conformal time rj. In the presence 
of small amplitude inhomogeneities the EoS can be decomposed into a homogeneous part and a small perturbation: 

wd = w D (rj) + 5w d (t), x l ) . (14) 

To obtain the matching conditions on this hypersurface, it is better for us to consider the general form of the perturbed 
metric 

ds 2 = a 2 {ii){{\ + 2A)dri 2 - 2B A dx l dr) - [(1 - 2?/')% + 2E^]dx i dx i } , (15) 

where commas denote derivatives with respect to spatial coordinates. Only two of the four variables A, B 1 ip, E are 
physical. Under the coordinate transformation 

i] — > fj = r) + £ and 

x l -> x 1 = x l + e l , (16) 

these metric perturbations transform as 

A A = A - - £ ' , B B = B + £° - f ' , 

V- -> $ = ip + H£,° , E^E = E-£, (17) 

and the perturbation of the EoS transforms as 

Swd — > Swd = Swd — w'jj^ . (18) 

We will use the temporal gauge to obtain the matching conditions. In this gauge, the matching hypersurface S 
coincides with fj = const, and the equation of this hypersurface 

w5_d(t/, x l ) = const. (19) 

implies 

frw D = 0. (20) 

Hence the time shift is 

e - ^ , (2D 

w D 

but ^ remains arbitrary. The induced 3-metric of this hypersurface and its extrinsic curvature are expressed as 

Qij = a 2 [(l-2Vi)%+2S ii ] , (22) 
Kij = ^ (H - HA - 4,') + ~ (E' - B) tij , (23) 

respectively. 

The matching conditions tell us that the induced metric and the extrinsic curvature should be continuous across 
the surface, i.e. that [qtj]± = and [-K"ij]± = 0. For the background, this requires that the scale factor a and the 
expansion rate T~L are continues. And for the perturbations, one obtains 

[i>}± = [E}± = (24) 
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and 

[HA + 4>'}± = [E -B]± = 0, (25) 

where the notation 

[<ip]± = 4>+ - -0-etc. (26) 

(the subscripts + and — indicating the values of the quantity on the two sides of the boundary) has been used. 
Making use of the gauge transformations (|17[) and (|21[) , wc obtain the matching conditions for the perturbations in 
an arbitrary gauge, 

w' D 

[UA + 4, 1 + {W -H 2 ) 5 -^] ± = , 

[#-iJ + ^] ± = 0. (27) 

Specifically, in the Conformal Newtonian gauge used in this paper (B = E = and $ = A, & = tp) these conditions 
become 



[*]± = o 
,8w D , 



w 



, 



D 



[H$ + *' + (H'-H 2 )^]± = . (28) 

D 

When $ = vp, i.e. in the absence of shear perturbations, and dividing the matter contributions into that of dark 
energy and that of regular cold matter, the Poisson equation ([5} becomes 

= -4arG{p D [6 D + (1 + w D )^9 D ) + Pm [S m + (1 + w m )^9 m ]} , (29) 

CI K K 

where the subscript m denotes matter. 

The first matching condition in (|28[1 means that the combination Sjj + (1 + wd)^6d should be also continuous. 
Because the matching hypersurface is characterized by wd = — lj one gets the following matching condition for the 
energy density perturbation of dark energy 

[5 D ]± = . (30) 
Now we turn to the physical meaning of the second condition in (|2"8"]l . After simple calculations one gets 



Sw D 1 8p D _ c 2 sD -w D 3H(l + w D ) 



_ _ WdSd ) = ^^[S D + ' -oj 0d] + l£ (31) 

w D w D pd w D k k z 

and at the matching hypersurface this becomes 

SWD c1d + 1 S d + ^9 d . (32) 



k 

Both Sd and w' D are continuous, and w' D must be non-zero in order to obtain crossing. Thus, the matching condition 
[Swd/w' d }± = implies that the momentum density perturbation of dark energy is also continuous, i.e. 

[0d]± = . (33) 

Eqs. (j30f and (|33|) coincide with the assumptions (fT3| used in Ref. Another way to see that (|30| and (|33|) arc 
valid we see that if these matching conditions are satisfied, then all of the matching conditions l|28p are satisfied. 
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Now, with the method discussed in this paper we can perform a numerical calculation to see how large the con- 
tribution of the dark energy perturbation can be. We modified and extended the CosmoMC code by implementing 
the dark energy perturbations discussed in this paper and take s = 10~ 5 , then fit the parameters of the dark energy 
EoS (wo,wi) to the current data from CMB observations including the 7-year WMAP t emp erature and polarization 
power spectra [H3|, and small-scale CMB measurements from BOOMERanG (mJ, CBI Hi, VSA H3 and ACBAR 
[54j . from the Union2 SNIa data set and from BAO[56]. In order to show the importance of the dark energy 
perturbation we have done the calculations separately for the two cases including and switching off the dark energy 
perturbations. In Figure [1] we plot our numerical results. One can see the obvious difference between the two cases 
given by the red solid line and the black dashed line. This is because the late time ISW effect differs significantly 
when dark energy perturbations are considered, and the ISW effects plays an important role on large angular scales 
for the CMB and the matter power spectra [57| . 



1 1 1 1 1 r 
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W 

FIG. 1: Two dimensional constraints on (wo,wi) from current observations of CMB + SNIa + BAO. The red solid and black 
dash lines represent the 2 a limits for the two cases with and without dark energy perturbations, respectively. The star 
represents the best ht value. 



IV. SUMMARY AND DISCUSSION 



In this paper we have revisited the dynamics of cosmological perturbations of dark energy and paid particular 
attention to the case when the EoS crosses the cosmological constant boundary. Single held or single fluid models, or 
scenarios based on a parameterized EoS of dark energy with a single component cannot cross wd = — 1 because the 
perturbations are singular and unstable at this point. The quintom model is able to cross this boundary naturally, 
however it requires more degrees of freedom, and lessons learned when studying the transfer of fluctuations through 
non-singular bounces makes us expect that, on scales smaller than the time duration of the transition phase, the final 
fluctuations will depend on the details of the model. This makes it hard to obtain a simple data fitting prescription. In 
particular, the more parameters are introduced, the more computing time is required for the numerical calculations. 
To obtain a simple way of analyzing data and assessing the observational evidence for or against the equation of state 
of dark energy crossing the cosmological constant divide it is thus very useful to have a prescription which does not 
introduce new parameters. 

In this paper we have presented a new approach to studying dark energy perturbations in the time interval when 
[ w + = —1 + e, W- = — I — e ] , i.e. during the crossing of the boundary wd = —1- We have proposed to apply the 
general relativistic matching conditions of [36l l37j . These conditions imply that the dark energy perturbations match 
continuously on the two sides of the surface wd = — 1- 

Let us mention some caveats to our analysis: Our method is applicable in the form presented here only if on either 
side of the matching surface all except for one fluid are negligible. Since at the crossing region this assumption will 
fail, this criterium implies that e cannot be too small. Secondly, since the dark energy fluctuations diverge when the 
equation of state crosses the cosmological constant line, then, in order to stay within the realm of applicability of 
linear cosmological perturbation theory, we have a second reason why e cannot be taken to be too small. On the other 
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hand, for length scales smaller than et, where t is the time when the EoS of dark energy crosses the cosmological 
constant line, the way in which the fluctuations pass through the transition region may depend on the specific quintom 
models. This argument prefers a small value for e. 

However, let us consider models where our assumptions are satisfied and where e is sufficiently small. Then, our 
results coincide with those of Ref. [H, and the arguments in this paper justify the method used in the numerical 
calculations of Since linear perturbation theory will break down as e — > 0, the small positive parameter e would 

not be taken to be too small. In Ref. [l[ it has been checked that with e ~ 10~ 5 , linear perturbation theory will be valid 
and at the same time the approximation of taking Sd and 6r> to be constant in the interval [w+ = — 1 + e, W- = — I — e 
] will yield results in agreement with those obtained by using actual quintom perturbations. 

Finally, with the method outlined in this paper we have performed a numerical determination of cosmological 
parameters. Our numerical results show explicitly the significance of the dark energy perturbations. 
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